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The NNLO QCD analysis of the CCFR data for xF^ : is there still the 
room for the twist-4 terms ? 
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The results of the more detailed NNLO QCD analysis of the CCFR data for xFa SF are presented. The 
factorization scale uncertainties are analyzed. The NNLO results for as{Mz) and twist-4 contributions are 
obtained. Despite the fact that the amplitude of the a;-shape of the twist-4 factor is consequently decreasing at 
the NLO and NNLO, our new QCD analysis seems to reveal the remaining twist-4 structure at the NNLO level. 
The definite N"^LO uncertainties are fixed using the [0/2] Fade resummation technique. 



It is known that the CCFR collaboration pro- 
vided not long ago rather precise experimental 
data for xF^ SF of i^N DIS and extracted the 
value of as{Mz) using the NLO DGLAP anal- 
ysis [Q]. In its process the twist-4 contributions 
were taken into account using the infrared renor- 
malon (IRR) model of Ref.[ || with its parameter, 
fixed as 1/2 of the originally proposed one. 

In the series of the subsequent papers [ ||]-[ D 
we concentrated on the attempts to fit the CCFR 
data at the NNLO level with the help of the Ja- 
cobi polynomial-Mellin moments version of the 
DGLAP method [ §- [||, based on the following 
equation: 



N„ 



N, 



x'^ii-xf V e'^^^{x)x (1) 



n=0 



where h{x) is the twist-4 contribution. 

We used the results of calculations of the 
NNLO corrections to the coefficient functions[||] 
and the analytical expressions for the NNLO cor- 
rections to the anomalous dimensions of the non- 
singlet moments with n — 2, 4, 6, 8, 10 [ 0, sup- 
plemented with the given in Ref. [ ^ n = 3, 5, 7, 9 
similar numbers, obtained using the smooth in- 
terpolation procedure of Ref. [ |ll[ . 

Using the fits with free Jacobi polynomial pa- 
rameters a, /3, we found that their values a w 0.7, 

« 3 corresponds to the minimum in the plane 
(a,/3). The form of h{x) was fixed (1) through 
the IRR model of Ref.[ |j with its coefficient A2 
considered as the free parameter and (2) as the 
function, modeled by free parameters hi = h(xi), 
where xi are the points of the experimental data 
bining. The QCD evolution of the moments has 
the following form 
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M„(Q2) [As{QI)J AD{n,Ql)C("){Ql) 

where A^iQ^) = a,(Q2)/(47r) is the WS- 
scheme expansion parameter, AD{n, Q^) = 1 + 
p{n)As{Q'^) +q{n)As{Q'^)'^ + ... comes from the 
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Table 1 

The Ql dependence of [MeV]. LO* means that in the LO-fits NLO as is used; NLO* {NNLO*) 
indicates that in the NLO (NNLO) fits NNLO {N^LO) as is used. [0/2] marks the resuhs of the N^LO 
expanded [0/2] Pade fits with as defined at the N^LO. 
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NNLO 


293 ±30 
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326±36 
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NNLO* 


284 ±28 


303±31 


308±32 


316±33 


316±33 


315±34 


[0/2] 


293 ±29 


323±32 


330±35 


335±37 


326±36 


319±35 



expansion of anomalous dimesnion, C^^\Q'^) — 
1 + C(i)(n)A,(Q2) + C(2)(n)A,(Q2)2 + the 

coefficient function of M„((5^) and M„((5o) is 
parametrized at the factorization scale Q^. In 
the case of / = 4 numbers of flavours the numer- 
ical values of p(n), C'^'^\n) and C(^^(n) are 
given in Ref. [ ^. We will use the expansion of 
as through the powers of 1 / ln{Q^ / Ajj^) in the 
LO, NLO, NNLO and N^LO, which contains the 
4- loop term of the QCD /3- function [[l^. 

Here we complete previous analysis of the 
CCFR data of Refs.[ ||, performed in the case 
of Qo = 5 GeV^, by varying Qo in the wide re- 
gion. The fits were done for the CCFR data, cut- 
ten at > 5 GeV^, without twist-4 effects, but 
with target mass corrections included. The re- 
sults for A.^^ are given in Table 1 for different Qq. 
The LO and NLO results are stable to the choice 
of Qq. The results of the LO* fits are higher 
than the LO ones, and from this level taking into 
account of other perturbative QCD effects are de- 
creasing the values of aIIL. The NNLO results 

IvI S 

are sensitive to the variation of Qo- The values of 
A^IL become stable for > 20 GeV^ only. The 

MS ^ — J 

same effect is manifesting itself for the results of 
the N^LO fits. This effects can be related to 
the peculiar behavior of the NNLO perturbative 
QCD expansion of M2. Using the numerical val- 
ues of p(2), g(2), C(i)(2) and C(2)(2), given in 
Ref. [ ^ , we obtain 

AD{2, Q^)C^^\Q^) = 1 - 0.132^(g2) (3) 
-46.155^(Q2)2 4. ... 



Thus, it is safer to start the evolution from Qq = 
20 GeV^, where the numerical value of the 
contribution in Eq.(3) is smaller. 

The results of our new fits of the CCFR data 
with twist-4 contributions, fixed through the IRR 
model of Ref.[ ||, are presented in Table 2. 

Table 2 

The results for A^-}— and the IRR coefficient A, 

MS 2 

[ Gel/2] f^Qjj^ the CCFR'97 data, cutten at > 
5 GeV"^ in different orders for Q^=20 GeV"^. 

A^ (McV) A^(HT) xVpoii^ 
LO 433±52 -0.33±0.06 82.8/86 

NLO 369±45 -0.12±0.06 81.8/86 

NNLO 326±35 -0.01±0.05 76.9/86 

N^LO 340±37 -0.04±0.05 77.2/86 

When the twist-4 parameters are not taken into 
account, the effects of the NNLO corrections are 
smaller than in the case of our previous analy- 
sis of Refs.[ ^, H (see Table 1). However, they 
are still sizable in the case when twist-4 contri- 
butions are fixed through the IRR model. They 
have the tendency to make the value of A2 com- 
parable with zero. As the result, the NNLO value 
of A.^1L is the same in the cases of both neglecting 
and retaining twist-4 terms. 

At Fig.l we present the extraction of the x- 
shape of the twist-4 terms from the LO, NLO, 
NNLO and expanded Pade fits with Qq = 
20 GeV^ for the unfixed twist-4 contribution. For 
A^IL we got: 331±162 MeV (LO level), 440± 183 
MeV (NLO level), 372± 133 MeV (NNLO level) 
and 371 ± 127 MeV ( expanded [0/2] Pade) 
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Figure 1. The x-shape of h{x) extracted from the 
fits of CCFR'97 for = 20 GeV^. 

One can see that taking into account of the higher 
order perturbative corrections is decreasing the 
ampHtude of the variation of h{x). This obser- 
vation is in agreement with the results of Refs.[ 
H, ^ , obtained for the case of Qq = 5 GeV'^ . How- 
ever, the change of factorization scale allows to 
detect the remaining twist-4 structure even at the 
NNLO. It is relatively stable to the application of 
the [0/2] Fade approximations method. 

When the twist-4 terms are fixed through the 
IRR model we obtain 



as{Mz)\NLO 

as{Mz)\NNLO 

and 

asiMz)\NLO 

o^s(.Mz)\nnlo 



0.120 ±0.002(siai) (4) 
±0.005(syst) ± 0.004(t/i.) 
0.118 ±0.002(siai) 
±0.005(syst) ± 0.003(t/i.) 



= 0.123 



+0.008 
-0.010 



[stat) 



(5) 



±0.005(syst) ± 0.004(t/i.) 

0.12lti]:0°^{stat) 

±0.005(syst)±0.003(t/i) 



when the twist-4 terms parameters h(xi) are free. 
The systematic errors are taken from the CCFR 
experimental analysis and the theoretical (th.) 
ambiguties are dominated by the uncertainty in 



the choice of the matching point in the NLO, 
NNLO and N^LO variants of the A/5-matching 
condition [ |l^ , derived following the lines of Ref. [ 
p^ . It was estimated by varying 6-quark thresh- 
old from Mf, = TOf, to Mi, = 6.5m;, [ |l^ and is of 
over ±0.002. 
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